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ABSTRACT 

In supercri tical systems the design inlet and 
outlet pressures are maintained above the thermo- 
dynamic critical pressure P c . Designers rely on 
this simple rule of thumb to circumvent problems 
associated with a subcritical pressure regime nested 
within the supercritical pressure system along with 
the uncertainties in heat transfer, fluid mechanics, 
and thermophysical property variations. The simple 
rule of thumb is adequate in many low-power designs 
but is inadequate for high-performance turbomachines 
and linear systems, where nested two-phase regions can 
exist. Examples for a free-jet expansion with back- 
pressure greater than P c and a rotor (bearing) with 
ambient pressure greater than P c illustrate the 
existence of subcritical pressure regimes nested 
within supercritical systems. 

NOMENCLATURE 

C clearance, cm 

D diameter, cm 

D c density at thermodynamic critical point, g/cm^ 
i index indicating inlet number 

e eccentricity 

G mass flux, g/s cm^ 

G r reduced mass flux, G/G* 


G* thermodynamic parameter, yp C p c /Z c , 
6010 g/s cm^ for nitrogen 

L spacer length, cm 

i aperture (orifice, Borda) length, cm 

P pressure, MPa 

P r reduced pressure, P/P c 


, Ohio 44106 

R radius 
Rg gas constant 
S entropy, J/g K 

T temperature, K 

T r reduced temperature, T/T c 
w mass flow rate, g/sec 

Z axial coordinate, cm 

Z c compressibility factor, P/pR g T 

e circumferential coordinate, rad 

p density, g/cm^ 

n viscosity 

to angular velocity 

Subscripts 

c thermodynamic critical point 

e exit 

in inlet 

s static 

0 stagnation 

INTRODUCTION 

Many current operational supercritical systems 
have been designed on a simple rule of thumb - 
maintain the design inlet and outlet pressures well 
above the thermodynamic critical pressure P c to 
circumvent two-phase flow and other associated 
uncertainties in heat transfer, fluid mechanics, and 
thermophysical property variations. The simple rule 


1 



of thumb is adequate in many low-power designs but, 
as will be demonstrated, it is inadequate for high- 
performance turbomachines and linear systems. 

Global thermodynamics permits local events that 
are considered unstable or unattainable when viewed in 
terms of total process dynamics. It becomes a matter 
of refining one's view of the black box containing the 
system. 1 For example, local "violations" of thermo- 
dynamics are termed metastabilities. Classic among 
these are the work of Skripov (1) and applied work 
such as fracture, boiling, and condensing (2,3). More 
recently the metastable character of water was 
investigated (4). 

Thus the concept of a two-phase region nested 
within a supercritical region would be considered not 
only unstable but unrealistic. However, there is no 
reason for such an exclusion as it depends on the 
thermohydromechanics that exist within components (the 
black box) operating at pressures above the thermo- 
dynamic critical point. Such conditions do not 
violate any thermodynamic hypotheses but are a source 
of incompatibility for the practicing engineer. 

This paper demonstrates that two-phase regions can 
be nested within a supercritical pressure system. 
Furthermore conditions for the existence of nested 
zones are prevalent for sharp-edge inlets and for 
eccentrically loaded high-performance turbomachines 
operating above the thermodynamic critical pressure. 

SIMPLE GEOMETRIES 

Borda-Inlet Tube 

Within a single-phase fluid, nested two-phase 
flows near sharp-edge inlets are commonplace and are 
in many cases represented as cavitation; however, the 
ambient pressure is subcritical. Consider the pres- 
sure profile of Fig. 1 (5-7). The inlet is sharp, and 
the L./D can be as high as 120 (6). The working fluid 
is nitrogen (P c = 3.417 MPa, T c = 126.3 K, 

Dc = 0.3105 g/crn^). The geometry is given in 
Fig. 1(a). The system was well instrumented and care 
was taken to maintain an adiabatic environment. For 
this configuration the inlet pressure was much greater 
than P c , but the pressure just 0.54 mm (0.020 in) 
from the sharp-edge inlet was less than saturation 
based on stagnation conditions^ P s (Sg). Downstream 
there is an immediate recovery toward P $ (Sq), or 
pseudocri tical^ P(Sg), but full recovery could occur 
anywhere between this point and the exit as illustrated 
by the three different profiles. These profiles were 
achieved by two methods: inlet stagnation condition 

control, and exit backpressure control. 

Inlet control . With controlled inlet stagnation con- 
ditions the pressure at the exit remained subcritical, 
but the pressure within the tube downstream of the 
inlet was supercritical for isotherms 127.8 K and 
130.5 K (Fig. 1(b)). This fully defines a two-phase 
flow nested between two zones at supercritical pres- 
sure. These results are somewhat open to question^ 
as Sq > S c . However, Sq < S c for run 1133 


Perhaps as a fracture surface that when viewed from 
close up appears with discontinuities everywhere yet 
appears continuous and smooth when viewed from afar. 
2 p s (Sq) is the static pressure as a function of 
stagnation entropy. 

^Pseudocritical : extrapolation of the saturation 

locus to the knee of the Joule-Tompson locus con- 
strained by the locii of specific heat maximums. 

^For Sq > S c , pseudonesting can occur, but two 
phases may not occur. 


(Fig. 1(c), where the recovery pressure was greater 
than Pp. Run 1133 is inlet temperature controlled 
and nested, but the nesting is marginal. 


Backpressure control . For run 1248 (Fig. 1(c)) the 
backpressure was increased sufficiently to move the 
pressure discontinuity ("shock") within the nozzle. A 
two-phase zone occurred but was at most weakly nested 
(all downstream pressures were less than P c ). For 
run 1547 the backpressure was increased above P c , but 
there was little change in the two-phase zone within 
the tube nested between Pq (at -27 cm) and P e 
(at 0 cm), where Pq and P e are greater than P c . 


Four Sequenti al Axially Aligned Borda or Orifice 
Inlets 

A more complex flow is established by a sequence 
of orifices or Borda apertures (8,9). The pressure 
profiles of Fig. 2(a) are for N Bordas in sequence 
with Tq approaching T c and backpressure depend- 
ent on downstream losses. The pressure signature here 
illustrates a nested set of two-phase zones within 
Borda 2, between the supercritical pressure zones of 
Borda 1 and Borda 3. With the inlet stagnation con- 
ditions constant the backpressure to Borda 4 was 
varied until P e > P c within Borda 2; the two-phase 
zones remain nested between Pq and P e , which are 
both greater than P c . Run 2879 (Fig. 2), illustrates 
that such zones can also exist within Borda 1. 

Because of an experimental facility limitation on Pq, 
nested regions could be approached in Borda 3 but not 
fully established; exclusion is therefore not war- 
ranted for either Borda 3 or Borda 4. 

For orifice flows, which may never choke, the 
pressures in the nested subcritical zones are often 
higher than saturation pressure. These zones shift to 
higher pressures as backpressure is increased and may 
not exhibit two-phase behavior (Fig. 2(b), orifice 2). 


TURBOMACHINE CONFIGURATIONS 

Normally seals and bearings for high-performance 
turbomachines operate at small eccentricities. How- 
ever, during abnormal events or even the limit-cycle 
type of event, dimensionless eccentricities can 
approach 1 (rubbing). 

Seal s 

Figure 3(a) illustrates the geometry of a high- 
performance turbomachine seal such as for the space 
shuttle main engine (SSME) (8,10). The working fluid 
is parahydrogen (P c = 1.2928 MPa, T c = 32.976 K, 

Dr = 0.0314 g/cnw). The pressure profiles 
(Figs. 3(b)-(d)) are for a fully eccentric, nonrotat- 
ing shaft. The inlet pressure is greater than P c , . 
as is the pressure at the inlet to the third step 
(Fig. 3(b)). However, immediately downstream of the 
step the pressure on one "side" of the seal (e = 0) 
remains above P c , while the pressure on the opposite 
"side" (e = n) is nearly P s . Thus the circumferen- 
tial pressure profile possesses all states between 
Pin > Pc 0=0 to approximately P s < P c at 
e = tt . Furthermore the system backpressure could be 
increased above P c , demonstrating yet another nested 
distribution - now in both e and Z (Fig. 3(c)). 

To corroborate these results, the tests were repeated 
( e . g . , Fig. 3(d)). Thus with eccentric seals pressure 
gradients of sufficient magnitude can exist to nest 
two-phase zones both axially and circumferentially. 


Bearings 

Unlike those for a seal the axial pressure 
gradients for a bearing are usually small, but the 
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circumferential gradients are usually much larger. 5 
With the establishment of nested zones in seal flows, 
such nested zones should occur naturally in bearings. 
The two-phase zone downstream of the minimum clearance 
in a bearing is well known, and the pressure and tem- 
perature on the stator have been mapped for a low- 
performance machine (12). It is also noted in (12) 
that the two-phase zone began to collapse as system 
ambient pressure was increased. However, calculations 
(13) for a similar high-performance bearing in fluid 
hydrogen at high eccentricity (Table 1) reveal a two- 
phase zone nested within a system where Pin > Pe > p c 
(Fig. 4). Nested two-phase zones in a bearing can 
lead to the onset of instabilities, contrary to the 
axial flow seal. The onset of the nested zone depends 
on the bearing parameter, 6 but onset is prevalent 
at lower ambient pressure (e.g., one-tenth those of 
Table 1), high eccentricity, and high rotating speed 
and for large geometries. This effect remains to be 
verified. If it is correct, this demonstrates yet 
another method of achieving nested two-phase zones, 
those not requiring a sharp geometrical configuration 
to "trigger" the event, spontaneous nucleation within 
the metastable field being necessary and sufficient. 

SUMMARY 

Zones of two-phase flows have been shown to be 
nested within a supercritical pressure envelope. 

These zones are established by controlling inlet 
stagnation conditions, exit backpressure, or combina- 
tions. The zones can be "triggered" by discontinuous 
changes in geometry or by significant pressure gradi- 
ents. For high-performance turbomachine seals and 
possibly bearings, nested two-phase zones can be found 
in the circumferential and axial coordinates. The 
phenomenon in bearings remains to be demonstrated 
experimental ly. 
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TABLE 1. BEARING PARAMETERS FOR ANALYTICAL SOLUTIONS 


Geometry: 



Diameter, mm ( in. ) . . 
Length, mm ( in. ) ... 
Clearance, mm (in.) . . 
Eccentricity 


.... 25.4 (1) 
.... 25.4 (1) 
. 0.0254 (0.001) 

Hydrodynamic : 



Inlet pressure, P-j n , MPa (psia) . . 
Exit pressure, P e , MPa (psia) . . 
Maximum pressure, P max , MPa (psia) 
Minimum pressure, Pm-j n , MPa^psia) 
Inlet temperature, T in , K (°R) . . 

. . 13.8 (2000) 

. . 13.45 (1950) 

. . 27.16 (3938) 
. . 0.079 (11.5) 
20 (36) 

Speed, rad/sec (rpm) . 



. 7850 (75 000) 




rl° DIFFUSER 


2.83(1.115) 


P 1 P 2 
dUJJ 

rl 0 DIFFUSER 

/ 

S- 




J3.65 . 
(1. 44) 

p-^?~ ^ 

- J '-0. 32 (0. 125) 

- — 0. 64 (0.25) 

J P 5 

Vf 

’3 ' p - 


14-L/D TEST SECTION 
USING BORDA INLET 


r 


0.08 

(0.03)1 



/Pi P 2 / 


STANDARD 3/4 AN 
FITTING 

^STANDARD 
■' 1/4 AN COUP- 
LING WITH Ni 
CONICAL 
WASHER 


MACHINED TO 
3/4 AN STANDARD 
FITTING^ 


ff 


S SL- 
TT“ 


"THROAT" 

0 

p 13 


-0.64(0.25) p ' 1 1 ' . 

— 0.65 (0. 255) P 4 6 8 10 16 P 18 
P 5 P 7 P 9 p 12 p 17 

L H 


STTT-tOflow 


\f ° 2 

\ L Z 27 (0.895) 
\ 

U Z83 (1. 115) 


— ^jp — 0. 04 - 0. 05 (0. 016 - 0. 02) 


0. 08 (0. 03) 

\ 

^-PRESSURE TUBE 
TO TRANSDUCER 

TYPICAL PRESSURE TUBE INSTALLATION 


(a) Configuration geometry. 

Figure 1. - Pressure profile for nitrogen flow through high-L/D tube 
with Borda inlet (Dimensions are in centimeters (inches). ) 
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RUN 

MASS 

INLET 

INLET 



FLOW 

PRESSURE, 

TEMPERATURE, 



RATE, 

p in» 

T in- 



w, 

g/sec 

MPa psia 

K 


1302 

1025 

10.54 1528 

130.5 

a 

1256 

1052 

10.45 1515 

127.8 

Q 

1294 

1122 

10.43 1512 

121.7 

□ 

1161 

821 

6.02 873 

116.8 

0 

1160 

893 

6.90 1000 

117.0 

A 

1159 

929 

7. 32 1 062 

116.0 

A. 

1249 
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86.2 


£ B 10. 45-MPa (1515-psia) STAGNATION ISOBAR 


4 A j 
0 O 

6 3-°l 


117 K STAGNATION 
ISOTHERM 


127.8 


r \0. 45-MPa (1515-psia) 
/// STAGNATION 
/ /' ISOBAR 


*— ^ o 


£ 10 



(b) Controlled inlet stagnation conditions. 
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(c) Controlled backpressure. 

Figure 1. - Concluded. 
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(a) Controlled inlet stagnation conditions, Borda inlets. 



(b) Controlled backpressure, orifice inlets (14). 

Figure 2. * Pressure profile for flow through four sequential axially aligned inlets. 
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(d) Backpressure control, repeated. 
Figure 3. - Concluded. 
















Figure 4. - Analytical pressure map for high-speed journal bearing 
with P j n > P e > P c . P jn * 13. 8 MPa (2000 psi); P e = 13.45MPa 
(1950 psi); N - 75 000 rpm ; = e v = 0. 68 ; Tj n * 20 K (36 °R) ; 

P b ' o. 1 MPa (!4. 7 psi) ; r i 56 K (100 °R). 
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